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ABSTRACT 

Supernova 2004dj in the nearby spiral galaxy NGC2403 was detected optically in July 2004. 
Peaking at a magnitude of 11.2, this is the brightest supernova detected for several years. Here 
we present Multi-Element Radio Linked Interferometer Network (MERLIN) observations of this 
source, made over a four month period, which give a position of R.A. = 07^37™17.044®, Dec = 
-|-65°35'57'.'84 (J2000.0). We also present a well-sampled 5 GHz light curve covering the period 
from 5 August to 2 December 2004. With the exception of the unusual and very close SN 1987A, 
these observations represent the first detailed radio light curve for the prompt emission from a 


Type H-P supernova. 

Subject headings: supernovae: individual (2004dj) 

1. Introduction 

We are currently using MERLIN^ (Multi- 
Element Radio Linked Interferometer Network) 
and NRAO’s VLA^ (Very Large Array) to mon¬ 
itor the radio emission associated with a sample 
of nearby starburst galaxies (Argo et al. 2004a) 
with the objective of determining their radio su¬ 
pernova (RSn) rates. As the active star-forming 
regions of these galaxies contain large numbers of 
massive stars, it is likely that most of the RSn 
will be associated with core-collapse supernovae 
(Type Ib/c and Type H). In addition to regular 
monitoring of the galaxies in our sample, we plan 
to use MERLIN to monitor the radio flux-density 
of other nearby supernovae events where possible 
(e.g. Beswick et al. 2004a,b) , particularly when 
the VLA is in compact configurations. 


^MERLIN is a national facility operated by the University 
of Manchester on behalf of PPARC in the UK. 

^The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under coopera¬ 
tive agreement by Associated Universities, Inc 


galaxies: starburst — radio continuum: supernovae 

The detection of SN 2004dj, visually on 31 
July 2004 (Nakano et al. 2004) and in the ra¬ 
dio using the VLA at 8.4 GHz on 2 August 
(Stockdale et al. 2004), provided an ideal oppor¬ 
tunity to study a nearby supernova and its evo¬ 
lution at radio wavelengths. The high declination 
of NGC 2403, the host galaxy (-1-65°), makes it an 
ideal source for high resolution MERLIN observa¬ 
tions. 

Peaking at magnitude 11.2, SN 2004dj was 
the brightest optical supernova for several years. 
Shortly after the initial detection it was reported 
that the spectrum of SN 2004dj showed features 
typical of Type H-P supernovae (Patat et al. 2004) 
Type H-P supernovae are believed to originate via 
core-collapse in hydrogen rich, massive stars. The 
main characteristic of Type H-P supernovae is 
that their optical luminosity, unlike other Type H 
supernovae, does not decay rapidly after peaking 
but remains approximately constant for two to 
three months before decaying exponentially. The 
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physics of supernovae are discussed in detail in 
Woosley & Weaver (1986) and references therein. 

Despite the fact that Type II-P supernovae are 
relatively common optically, few have been de¬ 
tected in the radio. As far as we are aware, only 
two Type II-P supernovae have been detected at 
radio wavelengths: SN 1987A (Ball et al. 1995; 
Kirk et al. 1994; Turtle et al. 1987, and refs, therein) 
and SN 1999em (Pooley et al. 2002). A further 
possibility is SN 1923A, although this is arguably 
a supernova remnant rather than a radio super¬ 
nova (Eck et al. 1998). The only detailed radio 
study of a Type II-P RSn is SN 1987A (see sec¬ 
tion 3.3), whereas the radio observations of SN 
1999em are rather limited with relatively few de¬ 
tections between day 34 and 70, despite extensive 
monitoring between 8.4 GHz and 1.4 GHz. Gon- 
sequently no well-sampled radio light curve of SN 
1999em was established. 

SN 2004dj occurred in NGC2403, a nearby 
(3.2 Mpc, Karachentsev et al. 2004) spiral galaxy 
located within the M81 group. As this is approx¬ 
imately half the distance of SN 1999em, this of¬ 
fered an excellent opportunity to establish a Type 
II-P radio light curve, so using a sub-set of the 
MERLIN array, we began monitoring this source 
frequently at 5GHz on August 5th and continued 
through to December 2nd 2004. 

2. Observations 

Monitoring observations at 4994 MHz with a 
bandwidth of 14.75 MHz were made with a sub¬ 
set of the MERLIN array (Thomasson 1986) be¬ 
tween 5th August and 6th October 2004. These 
employed a single 218 km baseline between the 
Cambridge and Defford antennas. Primary flux 
calibration was performed using the strong source 
3C84 assuming a flux density of 14.53 Jy, which 
was derived via a comparison with 3C286. The 
nearby point source 0733-1-646 was used for pri¬ 
mary phase calibration and for secondary flux cal¬ 
ibration assuming a flux density of 0.3566 Jy. Two 
additional secondary phase calibration sources 
(0713-1-669 and 0752-1-639) were also observed in 
order to establish the precision to which point 
source positions could be derived from single base¬ 
line data. Tests using all three phase calibration 
sources have shown that we are able to derive the 
position of SN 2004dj to better than 50 mas. 


The single-baseline ‘maps’, shown in Pig. 1, il¬ 
lustrate that SN 2004dj is detected at greater than 
3 ct in each of our 2.5 day integrations and that we 
are able to measure its position from such images. 
It should be noted that, although these single¬ 
baseline ‘maps’ can only provide limited structural 
information regarding the source, its position and 
flux density can be obtained. The light curve 
from August to early October (Pig. 1) is derived 
from 2.5 day vector averaged points after rotat¬ 
ing the data from the nominal pointing position 
to the measured position for SN 2004dj. Addi¬ 
tional points on the light curve (November on¬ 
ward) are derived from observations made using 
the full MERLIN array and were calibrated and 
imaged using standard methods. 

3. Results Discussion 
3.1. Position and progenitor 

From these MERLIN observations, the position 
of SN 2004dj is found to be R.A. = 07^37“I7!044, 
Dec. = -b65‘’35'57"84 (J2000.0), with an abso¬ 
lute error of <50 mas. No proper motion is de¬ 
tected between any of our epochs. This position 
is within 0.4 arcsec of the optical position re¬ 
ported by Nakano et al. (2004) and agrees with 
the Chandra X-ray position reported by Pooley & 
Lewin (2004). However, the VLA D-configuration 
position reported by Stockdale et al. (2004) 
(R.A. = 7‘'37“16!916, Dec = -b65°35'56"97) is 
~1.2 arcsec away from the MERLIN, Chandra 
and optical positions. The error in this radio 
position is almost certainly the result of confusion 
from the extended radio emission from NGC2403 
(Filippenko et al. 2004, S. D. Van Dyk, private 
communication). 

This position measurement puts the supernova 
near to n2403-2866 (R.A. = 07^37“16L93, Dec 
= -|-65°35’ 57'.'7; Larsen & Richter 1999) a clus¬ 
ter also known as Sandage 96, and the luminous 
blue variable (LBV) candidate 7-3909 (R.A. = 
07h37“16L93, Dec = -b65°35'57"6), another sug¬ 
gested progenitor (Weis et al. 2004). Optical ob¬ 
servations suggest that the object is extended and 
therefore more likely to be a cluster than a single 
LBV (Maiz-Apellaniz et al. 2004). 
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3.2. Radio light curve 

Using a single MERLIN baseline, the 5 GHz flux 
of SN 2004dj has been monitored on a regular basis 
from the 5th August through to the 6th October 
2004. Additional flux measurements using the full 
MERLIN array were made on several occasions in 
late-November and early-December. 

These radio observations, initiated just five 
days after the initial optical detection by K. Ita- 
gaki (Nakano et al. 2004), represent the best sam¬ 
pled early radio light curve of a Type II-P super¬ 
nova made to date, with the obvious exception of 
SN 1987A. At each epoch SN 2004dj was detected 
and imaged. The radio light curve for SN 2004dj, 
derived from the MERLIN observations, is shown 
in Fig. 1. 

On 5th August the prompt radio emission from 
SN 2004dj was detected with a flux density of 
I.SmJy at 5GHz (Argo et al. 2004b). Follow¬ 
ing this initial detection, the flux of SN 2004dj 
was observed to reduce during the second half 
of August before reaching an approximately con¬ 
stant flux value of ~0.9mJy. The flux of SN 
2004dj showed little significant variation from 
0.9±0.4mJy over the following month before de¬ 
caying slowly from this radio plateau level. By 
December (day ~125) the flux had decreased to 
the point where the supernova was undetectable 
above the noise (arms ~ O.lmJybm”^, December 
map, not presented here). It is not known when 
the 6 cm peak occurred, but it is possible that 
the earliest of our observations coincided with, 
or closely followed, the peak radio flux at 5 GHz; 
placing limits on both the timing and flux of the 
peak of the prompt radio emission. 

If it is assumed that our initial observations, 
in early August, do coincide with the 5 GHz peak 
in the prompt radio emission from SN 2004dj, this 
implies a peak flux of S 5 GHz = 1.9±0.1 mJy bm“^, 
equivalent to Lgcmpeak ~ 2.45 x 10^® WHz“^. At 
this peak radio luminosity, SN 2004dj is, along 
with the Type H-P supernova 1999em, one of the 
least luminous radio supernova detected (with the 
exception of the unusual LMC supernova 1987A, 
which was about two orders of magnitude less lu¬ 
minous). 

Based upon radio observations of several Type 
H supernova, Weiler et al. (2002) derive an empiri¬ 
cal relationship between peak luminosity and time 


delay between explosion and its peak such that, 

Lscmpeak 5.5t®3,!i X lO^® 

X (teem peak -to)^-^^°-"W Hz-1 (1) 

where teem peak — to is the number of days be¬ 
tween the supernova detonation and the peak of 
the prompt radio emission at 6 cm. Assuming the 
peak luminosity of SN 2004dj at 6 cm is 2.45 x 
10i®(±10%) WHz”! this implies teem peak — to = 
15illo days, placing the detonation date between 
July 11th and the first detection of SN 2004dj, July 
31st. This is broadly consistent with spectroscopic 
observations of SN2004dj, obtained on August 
3rd (Patat et al. 2004), which report SN 2004dj 
to have a spectrum similar to that of SN 1999em 
about three weeks after explosion (i.e. around July 
13th). Supernova 2004dj is a Type H-P super¬ 
nova and only the second of this Type for which 
a well sampled radio light curve exists. Thus, al¬ 
though these observations appear to be approxi¬ 
mately consistent with equation 1, more observa¬ 
tions of the prompt radio emission from Type H-P 
supernovae are required to make conclusions re¬ 
garding the whole of this class of supernovae, if 
their radio emission is to be used as a reliable lu¬ 
minosity distance indicator. 

3.3. Comparison with SN 1987A and SN 
1999em 

The only other Type H-P supernova to have 
been extensively monitored and imaged at ra¬ 
dio wavelengths in its early evolution is the su¬ 
pernova 1987A (Ball et al. 1995; Kirk et al. 1994; 
Turtle et al. 1987, and refs, therein). By virtue 
of the proximity of SN 1987A, situated in the 
LMC at a distance of ~50kpc, the prompt ra¬ 
dio emission from this supernova was detected a 
couple of days after its optical detection and has 
been monitored at various frequencies ever since 
(Turtle et al. 1987). 

The peak luminosity of the prompt radio 
emission from SN 1987A at GHz frequencies 
was 4xlO^®WHz-^ (Turtle et al. 1987), approx¬ 
imately two orders of magnitude lower than the 
peak luminosity at 5 GHz recorded for SN 2004dj 
(^ 2.6 X 10^®WHz-^). However, the 5GHz peak 
luminosity of SN 2004dj is comparable to that 
of Type H-P supernova 1999em which was esti- 
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mated to have a Lg cm peak ~ 2.2 x 10^®WHz“^ 

(Pooley et al. 2002). The X-ray luminosity of SN 
2004dj is approximately three times greater than 
that of SN 1999em (Pooley & Lewin 2004). 

Our estimate of tgcmpeak — to = ISillodays, 
and the spectroscopically based estimates of 
~30days of Patat et al. (2004), are compara¬ 
ble to the estimates for SN 1999em (~34days) 
(Pooley et al. 2002). In both of these cases, the 
time delay before the peak of the supernova’s GHz 
radio emission is an order of magnitude longer 
than that observed for the Type II-P supernova 
1987A. 

Along with the detection of prompt radio emis¬ 
sion from supernovae 1987A, 1999em and 2004dj, 
all three of these sources have also been de¬ 
tected at X-ray wavelengths (Masai et al. 1987; 

Dotani et al. 1987; Rodin et al. 1987; Pooley et al. 2002; 
Pooley & Lewin 2004). This is often true for 
core collapse supernovae (see figure 1 of Immler 
& Lewin 2003) and is unsurprising since both 
types of emission arise from the shocked circum- 
stellar medium. Consequently, it appears that an 
X-ray detection of a Type II-P supernova may be 
an indicator that it can also be detected at radio 
wavelengths, or vice versa. However, to date very 
few supernovae of this type have been detected in 
either of these wavebands. 

Interestingly, in the only other case where a 
Type H-P supernova has been extensively mon¬ 
itored at radio wavelengths, SN 1987A, the ra¬ 
dio flux was found to evolve and brighten af¬ 
ter ~1400 days (Manchester et al. 2002, and refs, 
therein), signifying the birth of a radio remnant. 

In fact, since this re-detection of SN 1987A, its 
radio emission has continued to increase at an ap¬ 
proximately monotonic rate and now far exceeds 
the radio brightness of its initial prompt radio 
emission phase. This increase in the observed ra¬ 
dio luminosity of SN 1987A occured as the ex¬ 
panding shock front from the supernova encoun¬ 
tered a large increase in the density of the circum- 
stellar medium. Coincidently, this second phase 
of emission was also observable at X-ray wave¬ 
lengths. Uniquely, SN 2004dj’s proximity may 
provide a second opportunity to study the tran¬ 
sition of a radio supernovae into a supernova rem¬ 
nant, although this will be dependent upon the 
density of the surrounding circumstellar medium, 
and thus the nature of the progenitor. Addition¬ 


ally if, in a few years time, SN 2004dj does signif¬ 
icantly brighten, it may be detectable and resolv¬ 
able using high sensitivity VLBI observations. As 
such it is essential that continued radio monitoring 
observations of this source are made. 

We would like to thank all the MERLIN staff 
for their assistance in making these observations 
possible. MKA acknowledges support from a 
PPARC studentship. RJB would like to acknowl¬ 
edge financial support by the European Commis¬ 
sion’s 13 Programme “RADIONET” under con¬ 
tract No.505818. 

REFERENCES 

Argo, M. K., Muxlow, T. W. B., Pedlar, A., 
Beswick, R. J. & Strong, M., 2004a, MNRAS, 
351, L66 

Argo, M. K., Muxlow, T. W. B., Beswick, R. J., 
Pedlar, A., & Marcaide, J. M., 2004b, lAU Cir¬ 
cular 8399, 3 

Ball, L., & Kirk, J. G., 1992, PASA, 10, 38 

Ball, L., Campbell-Wilson, D., Crawford, D. F. & 
Turtle, A. J. 1995, ApJ, 453, 864 

Beswick, R. J., Muxlow, T. W. B., Argo, M. K. & 
Pedlar, A. 2004, lAU Circular 8332, 2 

Beswick, R. J., Muxlow, T. W. B., Argo, M. K., 
Pedlar, A. & Marcaide, J. M., 2004, lAU Cir¬ 
cular 8435, 3 

Dotani, T., Hayashida, K., Inoue, H., Itoh, M. & 
Koyama, K. 1987, Nature, 330, 230 

Eck, C. R., Roberts, D. A., Cowan, J. J., & 
Branch, D., 1998, ApJ, 508, 664 

Filippenko, A. V., Li, W., Challis, P. & Van Dyk, 
S. D., 2004, lAU Circular 8391 

Immler, S. & Lewin, W. H. G. 2003 in Super¬ 
novae and Gamma-Ray Bursters. Edited by K. 
Weiler., Lecture Notes in Physics, vol. 598, 
p.91-111 

Karachentsev, Igor D., Karachentseva, Valentina 
E., Huchtmeier, Walter K. & Makarov, Dmitry 
1. 2004, AJ, 127, 203 

Kirk, J. G., Duffy, P. & Ball, L. 1994, ApJS, 90,807 


4 



Larsen, S. S., & Richtler, T., 1999, A&A, 345, 59 


Maiz-Apellaniz, J., Bond, H. E., Siegel, M. H., 
Lipkin, Y., Maoz, D., Ofek, E. O. & Poznanski, 
D., 2004, ApJ, 615, L113 

Manchester, R. N., Gaensler, B. M., Wheaton, 

V. C., Staveley-Smith, L., Tzioumis, A. K., 
Bizunok, N. S. Kesteven, M. J. & Reynolds, 
J. E., 2002, PASA, 19, 207 

Masai, K., Hayakawa, S., Itoh, H. & Nomoto, K, 
1987, Nature, 330, 235 

Nakano, S., Itagaki, K., Bouma, R. J., Lehky, M. 
& Hornoch, K, 2004, lAUC Circular 8377, 1 

Patat, F., Benetti, S., Pastorello, A. & Filippenko, 
A. V., 2004, lAU Circular 8378, 1 

Pooley, D., & Lewin, W. H. C., 2004, lAU Circular 
8390, 1 

Pooley, D., Lewin, W. H. C., Fox, D. W., Miller, 

J. M., Lacey, C. K., Van Dyk, S. D., Weiler, 

K. W., Sramek, R. A., Filippenko, A. V., 
Leonard, D. C., Immler, S., Chevalier, R. A., 
Fabian, A. C., Fransson, C. & Nomoto, K., 
2002, ApJ, 572, 932 

Rodin, V., Reppin, C., Pietsch, W., Engelhauser, 
J., Trumper, J.,Voges, W., Kendziorra, E., Be¬ 
zier, M., Staubert, R., Brinkman, A. C., Heise, 
J., Mels, W. A., Jager, R., Skinner, C. K., Al- 
Emam, O., Paterson, T. C. & Willmore, A. P. 
1987, Nature, 330, 227 

Stockdale, C. J., Sramek, R. A., Weiler, K. W., 
Van Dyk, S. D., Panagia, N., Pooley, D., Lewin, 

W. & Marcaide, J. M., 2004, lAU Circular 8379, 
1 

Thomasson, P., 1986, QJRAS, 27, 413 

Turtle, A. J., Campbell-Wilson, D., Bunton, J. D., 
Jauncey, D. L., & Kesteven, M. J., 1987, Na¬ 
ture, 327, 38 

Weiler, K. W., Pangia, N., Montes, M. J. & 
Sramek, R. A., 2002, ARA&A, 40, 387 

Weis, K., Romans, D. J. & Close, K., 2004, lAU 
Circular 8384, 4 

Woosley, S. E. & Weaver, T.A., 1986, ARA&A, 
24, 205 


This 2-column preprint was prepared with the AAS lAT^ 
macros v5.2. 


5 



1 1 1 

11th August 2l)04 

r 

n 1 1 

- <s 

4th September 2004 

, o 

r 

1 1 1 1 1 j 1 1 1 

23rd Novemb,er,2004 ® 

C; , ei‘ 

© - 


' 

0 






- '' 


. - © 
t> s 


? 0 
k) 


© 



^ '' Q ■ 





(5 Q 

^ 0 ^ 

0.54mJy/bm o 

1-'.1 1 1 1 

Map peak 1.21mJy/bm 

_1_Clj_1_ 

g] 

0.79mJy/bm ^ 

1 __3_1_ 

g 


o 

o 

58.2- 

o 

CM 

Z 

58.0 

o 

H 

< 

57.8 

Z 

-1 

O 

57.6 

LU 

a 

57.4 


07 37 17.10 17.05 17.00 

RIGHT ASCENSION (J2000) 


07 3717.10 17.05 17.00 07 37 17.10 17.04 17.00 

RIGHT ASCENSION (J2000) RIGHT ASCENSION (J2000) 



Day(0=31st July 2004) 


Fig. 1.— Radio light curve of SN 2004dj from MERLIN data. Day 0 of the radio light curve is 31st July 2004. 
The three upper panels show contoured images of SN 2004dj at different stages of its evolution. Observations 
up until October were performed using a single 217 km baseline, following this the entire MERLIN array was 
used. Each image is contoured with multiples of times 165, 63 and 132/iJy beam“^, from left to right 
respectively. The convolved beam size in the three images is 45x45 mas, 45x45 mas 58x57mas respectively. 
The peak flux of the source at each epoch is shown in the bottom left-hand corner of the images. The dates 
shown on each of the three images represent the start date for each of the imaged epochs. 
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